Collapse and Revival of 'Schrodinger Cat' States 
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We study the dynamics of the Jaynes-Cummings Model for an array of A^, two level systems (or 
qubits) interacting with a quantized single mode electromagnetic cavity (or quantum bus). For an 
initial cavity coherent state \a) and the qubit system in a specified 'basin of attraction' in its Hilbert 
space, we demonstrate the oscillation of a superposition of two macroscopic quantum states between 
the qubit system and the field mode. From the perspective of either the qubit or the field system, 
there is collapse and revival of a 'Schrodinger Cat' state. 



Quantum superposition of macroscopically different 
states of matter are of fundamental conceptual interest in 
many fields of physics such as Measurement Theory [l[ , 
Quantum Optics 0, Macroscopic Quantum Tunnelling 
[3| and Quantum Computation [J| . Many of the relevant 
experiments revolve around the preparation of and mea- 
surement on a specific class of states called 'Schrodinger 
Cat' states These are quantum superpositions of 

states which correspond to two (or more) different values 
of a macroscopic variable, such as the magnetization or 
the electric field in the cases of a large spin cluster, or 
large photon number or phase in a cavity mode. In this 
letter we report the surprising discovery that in a Jaynes- 
Cummings model (JCM), which describes an array of 
two-level systems (qubits) interacting with a single cav- 
ity mode, the time evolution can be such that sometimes 
the radiation field and sometimes the qubit subsystem is 
in a 'Schrodinger Cat' state. In other words, a superpo- 
sition of macroscopically different states can shift from 
one set of physical variables to another, as a function of 
time. 

The quantum dynamics of two- level systems (qubits), 
coupled to a single mode of an electromagnetic cav- 
ity, arise in many different physically interesting sys- 
tems. These include Rydberg atoms |5[, NMR studies of 
atomic nuclei [1,0], Cooper Pair Boxes Hi , Cavity Quan- 
tum Electrodynamics 0] , trapped ions [lo| and Quantum 
Computing [4|. A very general and simple Hamiltonian 
that captures the relevant physics in all these fields is the 
JCM [ll| (for one qubit) and its generalization for multi- 
qubit systems by Tavis and Cummings [l^ . Thus, our 
results are pertinent to a broad range of physical systems. 

One of the most interesting and surprising predictions 
of the JCM is the 'collapse and revival' of Rabi oscil- 
lations of the occupation probabilities for various qubit 
states as the system evolves, from an initial state which 
is a product of a coherent state \a) for the radiation field, 
and a generic qubit state jV'-A?,) HIj where Nq is the num- 
ber of qubits. These remarkable dynamics occur only 
because both the matter and the cavity field are treated 
fully quantum mechanically. Indeed, our aim here is to 
study the 'collapse and revival' of 'Schrodinger cat'-like 



states in a multi-qubit subsystem, as well as in the cavity 
field . 

For clarity let us specify the multi-qubit JCM Hamil- 
tonian, where each qubit labelled i has ground (excited) 
state \gi) (le^)) with energy eg^i (ce.i)- Up to a constant, 
the Hamiltonian has the form 
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where = |e,)(ei| - \gi){gi\, af = \ei){gi\, = \g^){e^\ 
and (d) is the creation (annihilation) operator for a 
photon with frequency w. The cavity-qubit^ coupling 
constant is \i and h£li = e^.i — eg^i. Here we consider 
only the cases of resonance, so w = fi^ for all i, and 
uniform coupling, so Ai = A. 

The celebrated 'collapse and revival' can be observed 
in the one-qubit case 0. It follows from an ini- 
tial system state of |5'i(0)) = |?/'i)|a), where |a) = 
g-|"l a = VWe^*^ , n is the average 

number of photons in the field and = {Cg\g)+Ce\e))- 
The Rabi oscillations (in the probability of the qubit be- 
ing in its initial state) demonstrate collapse, on a time 
scale oi tr ~ and then revival at tr ~ 221^^ xhis 



is illustrated in Fig. [T] for Cg = 1, Cg 
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by plot- 
ting Y^°^=o 1(5' '^l^i(^)) I I where {g, n\ corresponds to the 
qubit in its ground state with n photons in the cavity. 
A second notable feature of this time evolution, dis- 



131, is that at t 
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covered by Gea-Banacloche 

l'I'i(iir)) again factorises into a qubit part l^/')^^^ 
cavity part |$(itr)). Moreover, remarkably, the former 
is given by 
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where 9 is the phase of the initial coherent state, for all 



initial conditions such that I Co 



1 . Under 



same conditions |'!/')att attained sX t = ^tp- Because 
of this strikingly non-linear behavior, following Phoenix 
and Knight [IJ], we shall refer to the states Eq. ^ 
as 'attractors'. The probability that the qubit is in 
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FIG. 1: (color online) Time evolution for a system with one 
qubit. (a) the entropy of the qubit. (b) the probability of 
being in the qubit's initial state \g). (c) the probability of 
being in the state li/')^^^. At tr/2 the probability of being 
in the 'attractor' state goes to one while the entropy goes to 
zero. The qubit starts in the initial state \g) and the value of 
n = 50. 
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FIG. 2: (color online) Time evolution for a system with two 
qubits. (a) the entropy of the qubits. (b) the probability 
of the two qubit state \gg). (c) the probability of being in 
the two qubit 'attractor' state |'i/'2)^tt when the initial phase 
of the radiation field is 9 = 0. The two qubit 'attractor' 
state is reached at tr/4. The initial state of the qubits is 
-7=(|ee) + \gg)) and the value of n = 50. 



the State lij)^^^^ given by Y.n=o | (V'att, "-l^i W) | is also 
shown in Fig.[Tl together with the von Neumann entropy 
S"J(i) = -Tr(p9(<)lnp9(i)) associated with the qubit 
density matrix p'^(t) = Tr_F (|4'i(t))(*i(i)l), reduced by 
tracing over the field. Clearly, at t = the entropy 
S"'(i) approaches zero, and the analytical solution shows 
that in the limit n —^ oo the entropy goes to zero, in- 
dicating that the radiation field and the qubit are not 
entangled 

Prompted by these results, we have investigated the 
A^g > 1 qubit evolution, starting in a state |^'Ar,(0)) = 
I V'A", ) I Q^) ■ Analytically, in the large n limit, we have 
found that the states 
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can also be regarded as 'attractors' in a similar, dynam- 
ical, sense as outlined above. The only difference is that 
in the Nq > 1 case, I'^'Nq)^^^ only occurs at i = tr/2Nq 
for a restricted range of initial conditions, which we shall 
term the 'basin of attraction'. As before, at half way 
to revival of the initial state, namely t = tr/2Ng, the 
radiation field and the qubit system are not entangled. 
However, now the question of entanglement between the 
qubits arises. Compared to the one-qubit problem stud- 
ied by Gea-Banacloche [3] this is a new feature of the 
multi-qubit case. Clearly, as |i/'Af,)f(( is a simple product 
of individual qubit states when the qubit subsystem is in 
this state, the qubits are not entangled with each other. 
This is surprising because, as we shall show later, this 
attractor state can be reached from initial states with 
arbitrary entanglement between qubits. 

Before examining how such interesting dynamics can 
occur, it is useful to note that the above product state is 
a spin coherent state for a finite A^g-qubit system. Fol- 



lowing Radcliffe [l6[ we define such states as 
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where the states |7Vg,m) are the fully symmetrized Nq 
qubit states, for A^e qubits excited and Ng in the ground 
state, with m = and /3 a complex number 

that characterises the the state. The normalization 

and combinatoric factors are Af ^ [1 + ipr and 
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As can be readily shown. 



(^-m)!(^+™) 

the 'attractor' states identified in Eq. [3] are given by 
|/3 = ±ie**, iVg). Below we explore the implications of 
this observation for the dynamics of qubit states in the 
'basin of attraction'. 

Let us now investigate the 'basin of attraction' for the 
simplest multi-qubit system, with Nq = 2. In this case 
the time evolution described by 1^*2 (0) is readily found 
17[. For the most general, normalized, initial state 



IV2) = Cee \ee) + Ceg \eg) + Cq, \ge) + Cqg \gg) (5) 

the exact analytical solution will be given elsewhere [3 ■ 
Here we consider only the sector determined by the re- 
strictions: a = e'^Cee = ^~^^Cgg and \J \ ^ kl^ = 
Ceg = Cge- As wiU bc iUustratcd presently, these de- 
fine the 'basin of attraction' for the 'attractor' \'4'2)~^tf 
Namely, for any complex number a satisfying the condi- 
tion < |a| < l/x/2 in 



\^^)^a{e-'^'>\ee)+e^'\gg)) + 



\an\eg) + \ge)), 
(6) 
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FIG. 3: (color online) The qubit system started in the max- 
imally entangled state (|ee) + \gg))/y/2 and n = 50. (a) the 
entropy of the qubit system, (b) the probability of being in 
the state \gg). (c) the mixed state tangle of the qubit system. 



the probability (given by P2att{t) = {i^t att\ P'' {'t)\i^t att)) 
that the two qubits are in the state |'02)att approaches 
unity at some time t* . An example of such behavior (for 
6* = 0) is shown in Fig. [2l To highlight the similarity with 
the analogous phenomenon in the one qubit case (Fig.[T]), 
we also show the entropy S''{t). This is calculated from 
p'^, the two-qubit density matrix reduced with respect to 
the cavity field coordinate, which describes a mixed state 
for most times t. Notably, at t* = \tr, where P2att{t) = 
1, the entropy S''{t) approaches zero in the large n limit, 
indicating that the system of two qubits is not entangled 
with the field. 

As we have already noted, the interesting new feature 
of the two-qubit case as opposed to the one-qubit case 
is that the former is in general host to entanglement be- 
tween qubits and this provides an opportunity to study 
the dynamics of such entanglement. For example, whilst 
almost all of the initial states in the 'basin of attrac- 
tion' given in Eq. © describe entangled qubits, they all 
evolve into \4'2)au at t = jtj. where they are not entan- 
gled. We investigated the pure state tangle of the initial 
condition defined as r = 4 \CeeCg 



function of a and found that although there are only two 
points where r = 0, all values of entanglement, includ- 
ing T = 1 indicating maximal entanglement, are present 
in the 'basin of attraction'. Thus we are observing the 
time evolution of a generic amount of entanglement. To 
throw further light on the matter, in Fig. [51 w e show 
the time evolution of the mixed state tangle [20| calcu- 
lated from for the maximally entangled initial state 
\ip2) = (|ee) + \gg)). Evidently, just as the occupa- 
tion of the initial qubit states collapses and revives, so 
does the entanglement. This phenomenon was first noted 
by Rodrigues et al. in a similar context [21| . 

Interestingly, not only the 'attractor states', IV'a',)^^^, 
but also their 'basin of attraction' can be found for the 
general A'q-qubit case. In fact, using the large n expan- 



sion of Meunier et al. [2^ we found the 'basin of attrac- 
tion' for all the attractor states identified in Eq. [3l We 
have established that, for a given A'g, only initial states 
parametrized as follows: 



JV,/2 

E 

m=-Af,/2 



A(7V„a)e-^(^ 



if k is even 



-|iVg,77l) 



(7) 



if k is odd 



A{Nq,a) = 

where < |a| < , will evolve into the spin coher- 

ent state |/3 ~ ie^^^ ,Nq) at the time t — tr/2Nq. 

A remarkable feature of this 'basin of attraction' is 
that it includes 'Schrodinger cat' states of the qubit sys- 
tem. Indeed it can be readily shown, by carrying out the 
sum over m in Eq. ((7| and identifying two different spin 
coherent states, that 
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(8) 



Therefore the 'basin of attraction' consists of linear su- 
perpositions of two quite different spin-coherent states, 
Eq. H je-^^iV,) &nd\-e-''''Nq,). For Ng < 3 these 
two states are hard to distinguish, but for large Ng 
they can be regarded as macroscopically different and 
hence their superposition is very similar to the finite 
A^q-qubit 'Schrodinger cat' state with equal coefficients. 
Indeed for a = and a = ypj^, I^wJ^^q ^ 
^ mNg) + \-(3,Ng)) ^ \P,Nq),^^ for f3 = e^^ 

Thus, the dynamics governed by the JCM Hamilto- 
nian can transfer an initial state which is a product of a 
coherent state for the radiation field, \a), and a highly 
quantum mechanical 'Schrodinger cat '-like state for the 
qubits, at the time t* ~ iiito another product state 
where the qubit component is a rather classical qubit- 
(spin) coherent state, with no entanglement between 
the qubits. Given this intriguing time evolution, it is in- 
teresting to examine the state of the radiation field at t* . 
Again in the limit of large average photon number n in 
the cavity we find 
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V2^«-2 (9) 



Evidently, this is a 'Schrodinger cat '-like state for the ra- 
diation field. Hence, we may conclude that the initial 
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FIG. 4: (color online) Diagrams of the Q function [S^] (left) 
and spin Q function [2J| (right) at three different times, (a) 
the time t = Q, where the cavity is in a coherent state and 
the qubits are in a spin Schrodinger cat state, (b) the time 
t — tr/2Nq which is the time of the first attractor. (c) the 
time t = tr/Nq, when the field states are again overlapping 
and in a coherent state. They show the 'Schrodinger Cat' 
state moving from the qubits to the radiation field and back 
again in the limit n^oo, Nq =40, 6 = G. The Q function 
for the field has been scaled to a unit circle. 

quantum information which was encoded in the qubit 
'Schrodinger cat' state moves to and resides in the radi- 
ation field at t* . As might now be expected, the time 
evolution after t* returns this information to the qubits, 
as the initial state revives. This remarkable time evolu- 
tion is illustrated in Fig. [H 

In summary, we have studied the dynamics governed 
by the JCM for an array of qubits interacting with a 
single mode of quantized radiation field, a quantum bus. 
We have shown that a product state, comprising a coher- 
ent state for the photons and a macroscopic 'Schrodinger 
cat '-like state for the qubits, transforms into a product of 
a 'Schrodinger cat '-like state for the photons and a spin- 
coherent state for the qubits. Furthermore, as time goes 
on this transformation is reversed and then the process 
starts all over again. This suggests that some univer- 
sal quantum information is being passed back and forth 
between the two subsystems, each of which manifests it 
through its own physics. Further study of this intriguing 
'collapse and revival' phenomenon should include an in- 



vestigation of the effects of decoherence and dissipation, 
of which initial studies have been made by Mcunicr et al. 
[2^ . and design of experiments to observe it. 
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